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A class of cubic networks composed of a regular one-dimensional lattice and a set of long-range 
links is introduced. Networks parametrized by a positive integer k are constructed by starting from 
a one-dimensional lattice and iteratively connecting each site of degree 2 with a fcth neighboring 
site of degree 2. Specifying the way pairs of sites to be connected are selected, various random and 
regular networks are defined, all of which have a power-law edge-length distribution of the form 
P>{1) ~ with the marginal exponent s = 1. In all these networks, lengths of shortest paths grow 
as a power of the distance and random walk is super-diffusive. Applying a renormalization group 
method, the corresponding shortest-path dimensions and random-walk dimensions are calculated 
exactly for k = 1 networks and for k = 2 regular networks; in other cases, they are estimated by 
numerical methods. Although, s = 1 holds for all representatives of this class, the above quantities 
are found to depend on the details of the structure of networks controlled by k and other parameters. 

PACS numbers: 



I. INTRODUCTION 

Random walk in disordered environments is a much 
studied problem as it is a basic model in the theory of 
transport in heterogeneous media [7, and it has fea- 
tures much different from those characteristic of homoge- 
neous systems. In many cases, the diffusion law is altered 
in the way that the typical displacement of the random 
walker still grows as a power of time, i.e. R{t) ~ t", 
however, the diffusion exponent v differs from the value 
1/ = 1/2 characteristic of normal diffusion. In case of ran- 
dom walks with quenched random jump rates on regular 
lattices or random walks on fractal lattices (e.g. per- 
colation clusterslthe process is in general sub-diffusive, 
i.e. v < 1/2 P, Q]. The opposite case, > 1/2, where 
the particle is speeded up compared to normal diffusion, 
is much rarer. This phenomenon, called superdiffusion, 
arises e.g. in turbulent fluids [3], in chaotic transport 
in laminar fluid flows |^] or in systems of polymerlike 
breakable micelles Q. In random walk models, superdif- 
fusion can be induced by allowing long-range jumps: In 
the Levy flight 0, jumps of arbitrary length are executed 
with a probability that decays as a power of the length; in 
its "discretized" version, the Weierstrass walk, the pos- 
sible jump lengths are integer powers of an integer a > 1 

i- 

Another possibility for superdiffusion in case of random 
walks on a lattice is when long-range jumps occur only 
at certain sites of the lattice. In most cases, the under- 
lying lattice is a union of a regular and a random graph, 
where the latter has arbitrarily long edges. Such graphs 
have been investigated in various aspects. As a variant of 
the small-world network model constructed by rewiring 
edges of a regular network [3] , Newman and Watts have 
considered a one-dimensional lattice to which "shortcuts" 



'Electronic address: 'j uhasz@szfki.hu] 



between randomly chosen pairs of sites are added [sj . In 
this model, shortest pat hs M [§| , as well as spectral prop- 
erties of the Laplacian [l3| have been studied. Short- 
est paths or the diameter of d-dimensional lattices where 
long edges exist between pairs of sites with probabili- 
ties that decay algebraically with the distance have been 
studied b y se veral authors for d = 1 [ll|, [H, [H, [l3| and 
for d > 1 [l5|. Beside the above models, n-regular graphs 
of this type, where the degree of all nodes (i.e. the num- 
ber of edges emanating from a node) is n, have been in- 
troduced, as well. The issue of decentralized algorithms 
for finding short paths has been considered by Kleinberg 
in a d-dimensional lattice where each site has a fixed 
number of directed edges to randomly chosen sites that 
are selected with power-law decaying probabilities p^ . 
Benjamini and Hoffman have studied minimal paths in 
w-periodic graphs, which are unions of periodic graphs 
over the integers |T3] ■ Recently, Boettcher et al. have in- 
troduced a hierarchical 3-regular network consisting of a 
one-dimensional lattice and recursively constructed long- 
range links and calculated the random-walk dimension 
drw = l/i' of this network p^ . 

In this paper, we shall consider networks composed of 
a one-dimensional lattice and an additional set of long 
links. Particularly, we focus on the intriguing situation 
when the tail of the distribution of edge-lengths is of the 
form Py{l) ^ pi^-^. In this case, the diameter is conjec- 
tured to grow algebraically with the size of the network 
with a non-universal /3-dependent exponent p^ . [l5| . The 
aim of this work is to confirm this conjecture by explicit 
calculations and to probe whether the diffusion exponent 
in such networks displays a similar "marginal" behav- 
ior. For this purpose, we define a class of cubic (i.e. 
3-regular) networks with marginal edge-length distribu- 
tions (i.e. P>{1) ^ l^^) and shall constructively demon- 
strate that the shortest-path dimension that character- 
izes the size-dependence of the diameter and the diffusion 
exponent are not exclusively determined by the power in 
the edge-length distribution but depend on the details of 
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the structure of networks. The above two intrinsic prop- 
erties of networks are calculated exactly in certain cases 
by means of a renormalization group method; in other 
cases, they are estimated by the numerical implementa- 
tion of the renormalization procedure and by numerical 
simulation. 

The rest of the paper is organized as follows. In Section 
im a class of networks parametrized by a positive integer 
k is defined and its general features are discussed. In 
Section IIIIl networks with fc = 1 whereas in Section IIVI 
networks with fc > 2 are investigated in detail. Results 
are discussed in Section |V] and a heuristic derivation of 
the relation to resistor networks is presented in the Ap- 
pendix. 



II. GENERAL FEATURES 
A. Construction of the networks 

The networks to be studied in the subsequent sections 
have in common that they are all constructed in the fol- 
lowing way. A one-dimensional open or periodic lattice 
with N sites is given, where sites are numbered consecu- 
tively from 1 to N. The degree of all sites is thus initially 

2 in the periodic lattice and all but site 1 and N in the 
open one. Sites of degree 2 will be called in brief ac- 
tive sites. The links of this initial regular lattice will be 
termed short links in the foUowings in order to distin- 
guish them from long links generated by the following 
procedure. Let us assume that N is even and fc is a fixed 
positive integer. A pair of active sites is selected such 
that the number of active sites between them is fc — 1 
and this pair is then connected by a (long) link. That 
means, for fc = 1, neighboring active sites are connected, 
for fc = 2 next-to-neighboring ones, etc. This step, which 
renders two active sites to sites of degree 3, is then iter- 
ated until 2(fc — 1) active sites are left. These are then 
paired in an arbitrary way, which does not affect the ex- 
ponents appearing in asymptotic relations in the limit 
N oo. In the resulting network, all sites are of degree 

3 if the procedure starts from a periodic lattice whereas, 
in case of an open lattice, site 1 and site N remain of 
degree 1. The networks generated in this way (or, more 
precisely, the ensembles of networks in case of random 
networks) are characterized by the number fc and by the 
way pairs are selected. 

At some stage of the construction procedure, when the 
number of active sites is Na, these sites are distributed 
homogeneously on a coarse-grained scale ^ ^ 1/c, where 
c = Na/N is their number density. In other words, spac- 
ings between neighboring active sites have a rapidly de- 
caying distribution with the expected value 1/c. Thus, 
long links of length larger than I are generated typically 
when c is smaller than l/l and we obtain for the dis- 
tribution of lengths in an infinite network (N oo): 
Py{l) ^ l~^. Disregarding regular networks (see Section 
IIIIB II and IIVB[) , where exclusively long links of length 



In ^ (A > 1, n = 1, 2, . . . ) form, the probability of a 
long link of length I is thus inversely proportional to the 
square of the length: pi ~ l^^. Note, however, that not 
all edge lengths are realized in the construction proce- 
dure even for random networks. Namely, for fc = 1, only 
long links of odd length are produced. 

B. Studied quantities 

We are interested in two intrinsic properties of net- 
works. Beside the distance I measured on the underlying 
one-dimensional lattice, we also consider another met- 
ric: The chemical distance (or the length of the shortest 
path) £ between two sites is the minimum number of 
links that have to be traversed when going from one site 
to the other. We will see that, in the networks under 
study, the average length of shortest path between sites 
located in a dis tance I grows algebraically with I for large 
distances: ~ id-mia^ where the shortest-path dimen- 
sion c?min < 1 is characteristic of the particular network. 
This dimension describes at the same time the finite- 
size scaling of the diameter D{N) of a network, which is 
the maximum of chemical distances between any pairs of 
sites: D{N) - N'^--'-. 

The other quantity of interest is the random-walk di- 
mension of the network [l], [l^ . We consider a continu- 
ous time random walk on (infinite) networks, where the 
walker can jump with unit rate to any of the sites con- 
nected with the site it resides. The random- walk dimen- 
sion is defined through the relation [(a;^(i))]typ ^ 
^2/d,w^ where x{t) denotes the displacement of the walker 
(measured on the underlying one-dimensional lattice) at 
time t and [(a;^(i))]typ = expln(x^(i)) is the "typical 
value" of {x'^{t)). Here, (•) denotes the expected value 
for a fixed starting position in a fixed realization of the 
ensemble of networks, while the overbar stands for the 
average over starting positions and the ensemble of net- 
works. Note that the expected value {x^{t)) does not 
exist if t > since the expected value of edge lengths is 
infinite (in infinite networks). This accounts for that the 
average of ln{x'^{t)) is considered instead. The practical 
reason of the second averaging procedure is to eliminate 
random modulations which stem from the randomness of 
the structure of networks, whereas for regular networks, 
this averaging may be ignored. 

C. Relation to resistor networlcs 

The calculation of random-walk dimension is based on 
the well-known relation between the effective diffusion 
constant and the effective resistance of the equivalent re- 
sistor network [U, 01 • In the equivalent resistor network, 
each link has a (dimensionless) resistance related to 
the jump rate Wi along that link via — l/wi. Consid- 
ering networks built on an open lattice, the time t needed 
for the walker to get from one end of the network (site 
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1) to the other one (site N) is 

t - Nf (1) 

for large N, where f is the effective resistance of the 
equivalent resistor network between the two endpoints. 
For the sake of self-containedness, a heuristic derivation 
of this relation is given in the Appendix; for a precise 
formulation of this connection in an arbitrary network, 
the reader is referred to Ref. [l^l- Eq. ([1]) implies that 
the random-walk dimension is related to the resistance 
exponent C defined by the asymptotical relation r{N) ~ 
TV? as 

rfrw = 1 + C- (2) 

So, the problem of calculating dj-w is reduced to the cal- 
culation of the resistance exponent ^ of the equivalent 
resistor network. 

III. NETWORKS WITH k = 1 

The simplest class of networks is obtained when neigh- 
boring active sites are connected in the construction pro- 
cedure, i.e. k = 1. In this case, the elementary step of the 
calculation of the resistance exponent is that a minimal 
loop (pair of sites with a short and a long link between 
them) is eliminated as shown in Fig. [1] and the two sites 
next to the removed pair are connected by a single link 
with an effective resistance f. Applying the reduction 
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FIG. 1: Replacement of a minimal loop by a single link, 
rules of resistors in series and in parallel, we obtain 

r = Ta + Tc H . (3) 

rb + I 

Let us denote the chemical distance between site i and 
i + 1 hy ii. Initially, £i = 1 for all i and when the elim- 
ination step described above is carried out, the chemical 
distance transforms in a simple way: 

£ = ia + ec + i- (4) 

Now, a renormalization group scheme can be defined, 
in which minimal loops are eliminated one after the other 
as described above, exactly in the same order as the long 
links of loops were created in the construction procedure. 
As a consequence, the construction procedure and the 
renormalization can be performed simultaneously: once a 
pair of neighboring active sites is selected for getting con- 
nected it is immediately eliminated and replaced by an 
effective short link with resistance and chemical distance 
calculated according to Eq. ([3]) and Eq. ([5]), respectively. 
Let us start the construction-renormalization procedure 



with an infinite one-dimensional lattice {N oo) and 
assume that the length scale ^, which is the inverse of 
the number density c of active sites, is large, i.e. ^ ^ 1. 
Then the typical resistance and chemical distance on ef- 
fective short links scale with the length as r(^) ~ and 
£{0 ^ ^''™'», respectively. For large ^, the effective re- 
sistances and chemical distances are also large and the 
transformation rules read asymptotically as f ~ ra -I- Tc 
and i ~ £a + ic, where the relation / ~ g is meant as 
lim^^oo f / 9 — ^- As the transformation rules of r and £ 
are asymptotically identical, we conclude that for fc = 1 

C = rfmin. {k = 1) (5) 

Thus for fc = 1, both dj-w and dmin are determined by the 
resistance exponent C,. This exponent also has a further 
geometrical meaning for fc = 1 . Let us consider a network 
built on an open lattice and call a short link backbone 
link if its removal results in that the network becomes 
disconnected. Taking into account that shortest paths do 
not contain turnbacks (at least in the interior of the path) 
and they are composed of long links and backbone links 
alternately, one can see that the fraction of backbone 
links in a network of size N is proportional to N'^"^'"^ for 
large N . In other words, the set of backbone links is 
a fractal object characterized by the fractal dimension 

A. Uniform model 

Perhaps the simplest model in the class fc = 1 is ob- 
tained when pairs of neighboring active sites are selected 
equiprobably in the course of the construction procedure. 
We call this model the (fc — 1) uniform model. 

The resistance exponent of this network can be calcu- 
lated as follows. Consider an infinite system (TV oo) 
and assume that, at some stadium of the renormalization 
procedure, the number density of active sites is changed 
by an infinitesimal amount dc < 0. The differential of the 
"resistance density" p = re is then dp = cdf+rdc. On the 
other hand, dp ~ rdc/2, where the factor 1/2 comes from 
that in an elimination step two sites are deleted but the 
total sum of resistances is reduced only by r on average. 
Combining these equations, we obtain cdr ~ — r(ic/2, the 
integration of which results in the asymptotical relation 

Thus, in the fc = 1 uniform model, the dimensions under 
study are = dmin = 1/2 and d^w = 3/2. As can be seen, 
random walk is super-diffusive in this network. 

We mention that the calculation carried out above re- 
mains valid also in the general case when the jump rates 
are random variables provided the expected value of the 
inverse jump rates (i.e. resistances) exists. Thus we ob- 
tain the same random-walk dimension for such a disor- 
dered model, which shows that d^^ is determined solely 
by the structure of the network. 
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B. Closest- neighbor networks 

Another possibility to create networks with A: = 1 is 
when pairs are not selected with a uniform probabil- 
ity but always the pair (or pairs) of actually closest ac- 
tive sites are connected. These networks will be termed 
closest-neighbor networks. Here, distances between adja- 
cent sites are rendered initially unequal; they are either 
random or are modulated according to some aperiodic se- 
quence. Note that a periodic arrangement of short links 
of different lengths either would not produce arbitrarily 
long links or would result in equal distances between ac- 
tive sites at some stadium of the construction procedure, 
after which the procedure would no longer be unambigu- 
ous. 



1. Aperiodic networks 

First, we consider networks where the short links of 
different lengths are arranged according to aperiodic se- 
quences. These aperiodic sequences are composed of let- 
ters taken from a finite alphabet {a, 6, c, . . . } and are gen- 
erated by the repeated application of an inflation rule, 
which assigns a word (i.e. a finite sequence of letters) 
to each letter. For instance, the simplest sequence that 
is suitable for our purposes is the so-called silver-mean 
sequence. It is composed of two different letters, a and 
6, and is generated by the inflation rule a ^ Wa — aha, 
6 — > = a. Starting from letter a, the first few it- 
erations are a,aba,abaaaba,abaaabaabaabaaaba, etc. To 
these finite strings of letters, finite open lattices can be 
associated in which the two different edge lengths la and 
^h(< la) between adjacent sites follow the same sequence 
as the letters in the strings. In networks constructed 
in this way, there are multitudes of long links of equal 
length. These long links will be termed links of the same 
generation. In general, the aperiodic sequences that we 
need have to meet the following requirements. First, in 
order to avoid ambiguities in the construction procedure, 
adjacent short links of shortest length must not emerge. 
Second, the infinite system must be self-similar. That 
means when a new generation of long links is formed, the 
sequence of new distances has to follow the original aperi- 
odic sequence. Third, the order of distances must remain 
invariant when a new generation of links is formed, i.e. 
if la > lb holds, the new distances must satisfy la > h- 
If these requirements are fulfilled, a renormalization step 
in which a complete generation of minimal loops is elim- 
inated corresponds to a reversed inflation step. After 
the sporadic inventions of such sequences in the field of 
aperiodic quantum spin chains [21], an infinite class was 
introduced with the purpose of studying entanglement 
entropy in those systems [l^. In the inflation rule of 
these sequences, the words are composed of an odd num- 
ber of letters and letter b, which will represent shortest 
links by convention, stands at even places. The inflation 
rule of two-letter sequences with the above properties can 



be written in the general form 

( a^ Wa ^ aba{bay^^'^ , , 

"""'^ ■ \ ^ Wfc ^ a(5a)""\ ^ ' 

where n and m are integers fulfilling 1 < n < m. 
The special case m = n = I corresponds to the silver- 
mean sequence mentioned above while, with the choice 
m — n = 2, the well-known Fibonacci sequence is gener- 
ated. An example of three-letter sequences is the tripling 
sequence, generated by 

{a ^ Wa — aba 
b W}j = cbc (8) 
c ^ Wc = abc. 

Here, edge lengths are ordered as lb < Ic < la- The 
structure of a few aperiodic closest-neighbor networks is 
illustrated in Fig. O 




abaaaba abaabaa aba 




a babaabaa baba 




a bacbcabaa bccbcabca bacbcaba 



FIG. 2: Fragments of the silver-mean, the Fibonacci and the 
tripling network (from top to bottom). 

We shall show that the shortest-path dimension can be 
calculated exactly from the substitution matrix S of the 
underlying sequence. The elements Sap of this matrix 
are given by Sap = riaiuup), where naiwp) is the number 
of occurrences of letter a in the word wp. Consider- 
ing a string of letters Si and a column vector v(si) with 
the components [v{si)]a = na{si), it is easy to see that 
the application of the inflation transformation results in 
a longer string s^+i with the vector v{si^i) = Sv(si). 
Thus, the asymptotic ratio of lengths of successive strings 
obtained in the inflation procedure is given by the largest 
eigenvalue A+ of S. Besides, we need the asymptotical 
scaling of chemical distances between neighboring active 
sites when a new generation of links forms. The chemi- 
cal distance £a on an effective short link represented by 
letter a after a new generation of long links has formed 
is related to the previous chemical distances £a (o; 7^ b) 
as 

la =^^njiWa)i'y + nb{Wa)- (9) 

If the density of active sites is small, the chemical dis- 
tances are also large and the last term on the r.h.s. of 
Eq. ^ is negligible. We have thus the asymptotical rela- 
tion £a — X^-zs^fa '^7(^1)^7 ■ This yields that the chemical 
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distances grow in a renormalization step asymptotically 
by a factor A+, i.e. laj ia — ^^+1 ol ^ h, where A+ is the 
largest eigenvalue of the matrix S obtained from S by 
deleting the row and column related to letter h. As the 
length scale f grows by the factor in a renormalization 
step, we obtain finally that the shortest-path dimension 

is given by d,„i„ = 



and the random-walk dimension 



closely. In fact, the smallest value of djw that we found 
is drw = 1 + ~ 1.6309, which is reahzed in the tripling 
network. The geometry of this network is known to be 
extremal also with respect to the von Neumann entropy 
of aperiodic quantum spin chains 24] . In section IIVI we 



shall see that larger k may result in smaller random- walk 
dimensions. 



In Ah 



(10) 



For 
ing 



the family of networks constructed by us- 
two-letter sequences, we have A+(m, n) = 



4(m 



1) 



and A+(m, n) 



TO -|- 1. For example, for the silver- mean network we 



obtain d-c 



1 



1 



In 2 



In 3 



ln(l-|-y2) 

with the golden ratio 



for the Fibonacci network 
■ ^ In the 



tripling network constructed by using the rule in Eq. |[8]), 
the set of backbone links is closely related to the Can- 
tor set and accordingly, the random-walk dimension is 
= 1 + Numerical values of c?rw are shown in 

Table [T] for a few two- letter networks. It follows from Eq. 



m=n=l (silver mean) 


1.7864... 


m=n=2 (Fibonacci) 


1.7610... 


m=n=3 


1.7623... 


m=n=4 


1.7683... 


m=n=5 


1.7748... 



TABLE I: Random- walk dimension of networks constructed 
on the basis of the two- letter inflation rule in Eq. ([7]). 



pO)) that 1 < < 2 for the aperiodic networks intro- 
duced in this section. This means that random walk is 
super-diffusive in these networks just as in the uniform 
model. Next, we discuss the bounds of d^^ in this class 
of networks. As can be seen from the data in Tabled for 
the two-letter networks with n = m, the random- walk 
dimension increases monotonously with to (i.e. with the 
length of words) after the minimum at to = 2 and one 
can see from Eq. pU]) that it tends to 2 in the limit 
TO = n — *■ 00. This tendency is intuitively easy to under- 
stand: A fragment of these networks which corresponds 
to a single word has periodically arranged short links, 
where the links of shortest length 4 are located at even 
places. Therefore, the long links generated in the interior 
of such a fragment are of limited length (Zf,) and, as a con- 
sequence, diffusion in such fragments is normal. Thus the 
upper bound of drw in this family of networks is drw = 2 
and this value, which is characteristic of normal diffu- 
sion, can be approached arbitrarily closely by choosing 
sufficiently long words in the inflation rule. Concerning 
the lower bound of drw, presumably there does not exist 
representatives of the family of fc = 1 aperiodic networks 
which approach the ballistic limit dj-^ = 1 arbitrarily 



2. Random closest-neighbor network and related networks 

Another example for closest-neighbor networks is the 
one generated with random initial lengths. In the case 
when the initial distribution of lengths is discrete and 
there may be adjacent short links of shortest length 
present with finite probability, the construction proce- 
dure is extended with the additional rule that a pair is 
randomly selected from those with shortest distance with 
uniform probability. 

The renormalization group scheme of this network with 
the asymptotical rule £ ~ €a -|- £c is formally identical 
to that arising in the context of a model of coarsening 
introduced in Ref. 23]. For that recursion scheme, it 
has been shown that the exponent z that appears in the 
relation £(^) ^ between the variable £ and the length 
scale f is the zero of a transcendental equation and it 
has been found that z = 0.82492412 • • • . Thus for the 
random closest-neighbor network, we have d^^^ = z and 
dj-w = 1 + z. 

Next, we discuss a variant of the random-closest neigh- 
bor network, the renormalization of which is identical to 
that of certain quantum spin chains [24| and due to this 
equivalence, the random-walk dimension can be exactly 
determined again. Let us consider a one-dimensional lat- 
tice of size N with three variables at each short link: 
the length /„, the chemical distance In, both are initially 
equal to 1, and an independent, identically distributed 
random variable /?„ that we call /3-distance. Now, a cu- 
bic network is generated as follows. The pair of active 
sites with the shortest /3-distance, say /3(,, is chosen and 
the two sites are connected by a long link. The length 
and the chemical distance on the new effective short link 
produced in the equivalent renormalization step are cal- 
culated ordinarily but the /3-distance has an anomalous 
transformation rule: 



(11) 



where the indices refer to links as given in Fig. [TJ As 
the structure of the network constructed in this way is 
identical to that of singlet bonds in the so-called random- 
singlet phase of antiferromagnetic quantum spin chains, 
we call this network random-singlet network. In Ref. [2^ 1 , 
it was shown in the limit N —> 00 that, when the density 
c of active sites goes to zero, the typical value of the 
variable £ scales with the length ^ = 1/c asymptotically 

as £{^) ~ where = is the golden ratio. The 

shortest-path dimension of the random-singlet network 
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is thus dmin = 

dimension is dr 



w 0.8090, 

_ 5+\/5 
■ — 4 • 



and the random-walk 



Finally, we examine another variant of the random 
closest-neighbor network: This is constructed by con- 
necting pairs of active sites with the actually shortest 
chemical distance. This means physically that when 
searching the closest neighbor of an active site, the al- 
ready existing long links are also made use of. When 
generating this network, which we call random minimal- 
chemical-distance network, the initial chemical distances 
are made random; in the numerical calculations, we used 
the initial values In = 1 + e„, where £„ <C 1 is a small 
random variable. One can easily see that the variables 
In and in become positively correlated in the renormal- 
ization process of this network, i.e. for small values of In 
the variable ^„ is also typically small. On the grounds of 
this observation, we expect that the shortest-path dimen- 
sion of the random minimal-chemical-distance network is 
close to that of the random closest-neighbor network. Ac- 
cording to results of numerical calculations, this is indeed 
the case. We have calculated the effective resistance be- 
tween endpoints of networks generated from open lattices 
of size N — 2", n = A, ... ,14. Data obtained in this way 
in 10^ — 10^ independently generated networks were av- 
eraged for each N. Results are shown in Fig. [31 The 
method was tested on the uniform model, on the random 
closest-neighbor network and on the random-singlet net- 
work, for which we obtained C = 0.500(1), C = 0.825(1) 
and C = 0.809(1), respectively. The resistance expo- 
nent of the random minimal-chemical-distance network 
obtained in this way is C = 0.826(1), which is very close 
to that of the random closest-neighbor network. 
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FIG. 3: Effective resistance exponent C(-^) = 
ln[r(2A'")/r(A'^)]/ in 2 plotted against the system size for 
the = 1 uniform (□), the random closest-neighbor 
(o), the random-singlet (A) and the random minimal- 
chemical-distance network (v). Horizontal lines indicate the 
extrapolated values in the limit A'^ oo. Inset: Dependence 
of the average resistance on the size A*'. 



IV. NETWORKS WITH k > 1 

In this section, we will discuss networks with fc > 1, 
which differ from the class of networks with k — 1 studied 
so far in that the resistance exponent and the shortest- 
path dimension are no longer equal. We will focus mainly 
on the case k = 2, where a renormalization group scheme 
similar to that applied for k = \ can be formulated. 
Herein, the minimal loops, which are triangles for fc = 2, 
are eliminated one after the other as shown in Fig. [H The 




FIG. 4: Renormalization scheme for networks with k — 2. 

result of such a renormalization step is that site h and d 
together with the long link connecting them vanish. The 
new effective resistances can be calculated by employ- 
ing the star-triangle transformation of resistor networks. 
This yields 



nrf 


n 


+ rc 4 










n 








nrc 
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\-rc + 





(12) 



The new feature here compared to the renormalization of 
fc = 1 networks is that resistances of long links are also 
transformed. 

Note that a similar substitution by which the long links 
of minimal loops are removed without changing the topol- 
ogy of the rest of the network cannot be formulated for 
fc > 2. 



A. The uniform model 

First, the case is considered when pairs of active sites 
are selected randomly with uniform probability in the 
construction procedure. 

For fc = 2, when a triangle is eliminated in the renor- 
malization procedure, the change in the total sum of re- 
sistances (including those on long links) A can be written 
in the form 



-n 



7 



rb+rc + rf 



(13) 
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As can be seen, the average value of A cannot be ex- 
pressed by the average values of resistances r^jr^ and rj, 
therefore the evolution of the entire distribution of the 
latter quantities should be taken into consideration. An- 
other difhculty is that the renormalization rules in Eq. 
()12|) apparently induce correlations between the quanti- 
ties fa, fc and fg. Therefore, we resorted here to numer- 
ical methods again in order to estimate drw and dmin- 
According to results of numerical renormalization, the 
absolute value of the average of the last term on the r.h.s. 
of Eq. is growing proportionally to the average of 

effective resistances in the course of the procedure. Al- 
though, it is thus not negligible compared to the average 
of first two terms on the r.h.s. of Eq. (fT3|) . it is by 
an order of magnitude smaller than those two terms. If 
the last term on the r.h.s. of Eq. (|13p is omitted, the 
resulting simple renormalization process that is analyti- 
cally tractable may provide a first approximation for the 
resistance exponent. Using that the average resistance of 
short links is equal to that of long links, which follows 
after all from the symmetry of the transformation rules 
in Eq. (|12p. we can write for the differential of the re- 
sistance density in this simplified process: dp = irdc/i, 
which leads finally to r(^) ~ ■ 

We have carried out the renormalization procedure 
numerically until two active sites were left in networks 
generated from periodic chains of size N = 2'^, n — 
4, . . . , 15. The average resistance of the remaining two 
effective short links was calculated from data obtained in 
10^ — 6 X 10^ independent networks for each system size 
N. The resistance exponent extracted from the finite- 
size scaling of the average resistance is C = 0.311(1), see 
Fig. El 
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FIG. 5: Average resistance of the k = 2 uniform model plotted 
against th e system size . Inset: Effective resistance exponent 
({N) = ln[r(27V)/r(7V)]/ln2 plotted against the system size. 
The horizontal line indicates the extrapolated value in the 
limit N ^ oo. 

We have also performed numerical simulations of the 
random walk in k = 1,2,3 uniform networks of size 
N = 10^ and measured the displacement x{t) of the 
walker at time t = 2", n = 4, . . . , 15. As we measured 



x{t) in a given network iV-times using all the N sites of 
the network as starting positions of the walker, the aver- 
aging over different stochastic histories for a given s tart- 
ing position was ignored and the quantity ln|a;(t)| was 
calculated. In addition to the averaging over starting po- 
sitions, data obtained in 500 independent networks were 
averaged. Results are shown in Fig. [51 The estimated 
random- walk dimensions are d,-^(k = 1) = 1.49(1), 
d,^{k = 2) = 1.31(1) and d,^{k = 3) = 1.19(1). As can 
be seen, these values are less accurate compared to the 
resistance exponent calculated by numerical renormaliza- 
tion. Nevertheless, drw{k = 1) measured in this way is 
compatible with the exact value 3/2 and d^wik = 2) is 
compatible with ( calculated numerically for k = 2. 




FIG. 6: Time- dep endence of the average logarithmic dis- 
placement In \x{t)\ of the random walker measured in numer- 
ical simulations in the uniform model with different values 
of k. Inset: Effective random walk dimension drw(i) calcu- 
lated from neighboring pairs of data points as l/dr„(t) = 
[ln|a;(2f)| — In |3;(t)|]/ In 2. Horizontal lines indicate the ex- 
trapolated values in the limit N oo. 



Beside the random-walk dimension, we have measured 
the shortest-path dimension for k = 1,2,3, as well. The 
length of shortest-path between sites in a distance N/2 
was determined by a simple breadth-first search algo- 
rithm. This was performed in 5 x 10^ — 5 x 10'^ in- 
dependently generated networks for each size = 2", 
n = 4, . . . , 17, and for A^/8 pairs of sites in each network. 
The average chemical distance £{N), plotted against N 
in Fig. [71 was found to grow for large N a.s £{N) ~ 
TYdmin _|_ (J -vvith a constant term C « — 2, that shifts the 
effective dmin(A^) considerably for moderate N. Using 
this ansatz, the estimated shortest-path dimensions are: 
rfmin(fc = 1) = 0.500(1), d„,in(fc = 2) = 0.440(1) and 

rfmin(fc = 3) = 0.381(1). 
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FIG. 7: Average length of shortest paths 1{N) between 
sites in a distance I = N/2 plotted against A'^ for differ- 
ent values of k. Inset: Effective shortest-path dimension 
dniin{N) calculated from neighboring pairs of data points as 
dmin{N) = \n\I{2N/I(N)]/ln2. Horizontal lines indicate the 
extrapolated values in the limit A'' oo. 



B. Regular networks 

In addition to the random networks studied in the pre- 
vious section, one can define regular networks with fc > 1, 
as well. The cubic, hierarchical network studied in Ref. 
[isl l is an example for a regular network with fc = 2. 
In this section, we shall consider regular k = 2 networks 
which are defined by means of aperiodic sequences. Here, 
the random- walk dimension and the shortest-path dimen- 
sion can be calculated exactly. 

Let us consider the subclass of aperiodic sequences dis- 
cussed in section IIII B 11 where the length of words is 
either one or three. One can define networks by using 
inflation rules with longer words, as well, but we shall fo- 
cus on this simple subclass. Once a sequence of this class 
is given, a finite network with odd sites can be defined as 
follows. A finite string of letters which is generated from 
a single letter is taken and this time, not the links but 
the sites of a one-dimensional lattice are labeled with the 
letters of the string. The sites are grouped into blocks 
corresponding to words Wa in the inflation rule, which 
can be done unambiguously. Then, sites belonging to 
one-letter blocks are renamed according to the reversed 
inflation rule Wa a, where Wa is the one-letter word 
corresponding to the block. In blocks composed of three 
sites, the two lateral sites are connected, and the middle 
one is renamed again according to the reversed inflation 
rule Wa ct, where Wa is the word corresponding to the 
block. The above step is then iterated until only one ac- 
tive site is left. In finite networks this site is of degree 
2. 

We have not found a general way of calculating d^w 
and dmin in these networks so they have to be treated 
individually. Two examples will be discussed. 



1. The k = 2 tripling network 

First, we examine the k — 2 tripling network which 
is constructed by means of the tripling sequence with 
the inflation rule in Eq. ((S]). This network is shown in 
Fig. [H Let us generate a sequence of finite strings by 




FIG. 8: A finite k = 2 tripling network (top) and a finite 
k — 2 silver-mean network (bottom). 

applying the infiation rule n-times on letter a and con- 
sider the corresponding sequence of finite networks built 
on open lattices the sites of which are labeled by these 
strings. The size of the nth such network is Nn = 3". 
By constructing the first few such networks and observ- 
ing their self-similar structure, one can easily check that, 
for the length of the shortest path between the two end- 
points, the recursive relation in — 2£n-i + 1 holds. Thus 
for large n, in ~ 2". Expressing n with Nn yields 
£ ~ jYin2/in3^ rpj^^ shortcst-path dimension of this net- 
work is thus dmin = li:! ~ 0.6309. 

The resistance exponent of this network can be calcu- 
lated by a renormalization procedure in which the gen- 
eration of long links belonging to minimal loops is itera- 
tively eliminated according to the scheme shown in Fig. 
m Let us consider the infinite network {N —f oo), with 
resistance r on short links and resistance p on long links. 
When the generation of long links belonging to minimal 
loops is eliminated, we obtain a similar network but with 
modified parameters f and p. Using Eq. (jl2p , we obtain 
for the renormalized parameters: 



2rp 
2r +p'' 



p^p + 



2r2 



2r+p 



(14) 



These equations yield for the scaling factor of resistances: 
X = r /r = p/p = 5/3. Keeping in mind that the scaling 
factor of the length is A = 3, we obtain for the resistance 
exponent: C = = 
dimension: drw = 1 



_ In 5/3 £ random-walk 

In 3 
■5 In 3 



1.4650. 



2. The k — 2 silver-mean network 

Next, we discuss the k = 2 silver-mean network con- 
structed by using the inflation rule in Eq. ([7]) with 
m = n = 1 (see Fig. Observing the self-similar struc- 
ture of finite representatives of this network, one can for- 
mulate the following recursion equation for the lengths of 
shortest-paths between endpoints in finite networks con- 
structed by consecutive strings of letters: £k = 3£k-2 + 2. 
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The asymptotic scaling factor of the chemical distance is 
thus -x/S- As the asymptotic scaling factor of the size of 
the network is 1 + \/2 (see Sec. IIIIB II) . the shortest-path 



solves the above system of equations with the scaling 

1 



dimension is dm in = 



_ In \/3 



0.6232. 



ln(l + \/2) 

The calculation of for the silver-mean network is 
somewhat more involved than for the tripling network. 
In order to find blocks which transform in a self-similar 
way under a renormalization step, an additional site of 
degree 2 is inserted in the middle of each short link. The 
number of short links (and sites) is thus doubled but 
if we assign a resistance 1/2 to short links instead of 1 
in this extended network, the effective resistance is obvi- 
ously equal to that of the original network. Furthermore, 
instead of keeping track of resistances of long links, we di- 
vide the resistance of a long link equally among the two 
sites that it connects. This means formally that there 
is an additional variable at all sites of the original 
network (initially 1/2), which changes in the course of 
the renormalization procedure, just as the resistances of 
short links. After eliminating the first generation of long 
links, the effective resistances of short links emanating 
from sites of type a will be different from those of short 
links emanating from sites of type h. These resistances 
will be denoted by and r;,, respectively. The addi- 
tional variables will also be different at the two types of 
sites. These will be denoted by Pa and p;,. This structure 
of the parameters remains, however, invariant when the 
subsequent generations of long links are eliminated. The 
substitution of the two different types of building blocks 
of the silver-mean network is illustrated in Fig. [HI When 



Pa 



Pb 






FIG. 9: Renormalization scheme for the fc = 2 silver-mean 
network. Sites of the original network are symbolized by black 
circles while the additional sites by grey circles. 

in an arbitrary stadium of the renormalization procedure 
the parameters are ro,rf,,Po and then, after the next 
generation of long links is eliminated, the renormalized 
parameters can be expressed in terms of the original ones 
as 



Pa = Pb + 



Pa+ra+n 

{vq + nf 

'2{Pa +ra+ n) ' 



Pb = Pa- 



(15) 



After a lengthy but elementary calculation one obtains 
that the ansatz fa/va = h/n = Pal Pa = Pb/Pb = A 



factor A 



1+V2 



\/ 2y/2. The resistance exponent is 



thus C = j^Tj = In i^j^ + \/2V2j I ln(l -h V2) and the 
random-walk dimension of the fc = 2 s ilver-mean net- 
work is drw = 1 + C = In (l + V4 + 2V2j / ln(l + \/2) w 
1.4575. 



V. DISCUSSION 

We have seen that, in the networks studied in this 
work, the distribution of edge lengths is broad and, as 
a consequence, the length of shortest paths grows sub- 
linearly with the distance and random walk becomes 
super-diffusive. A model similar to random networks of 
this class (although, not n-regular) is the model by Ben- 
jamini and Berger, where any pairs of sites in a distance I 
are connected with probability "pi ~ f3l~^~^ for large I. In 
that model, the diameter D, which is measured in terms 
of chemical distances, was conjectured to grow with the 
system size as D{N) Ar'imi„(/5) ^itj^ 

some /3-dependent 

exponent < (imin(/3) < 1 in the marginal case s = 1 
|12 |. Later, a power-law lower bound for /? < 1 and a 
power-law upper bound for any f3 were shown to exist 
asymptotically with high probability for the size depen- 
dence of the diameter [1^. Thus, as far as the length 
of shortest paths (or the diameter) is concerned, random 
networks studied in this paper behave as that model is 
conjectured to do. Furthermore, we have shown that be- 
side dmin, the random- walk dimension is also influenced 
by the details of the structure of networks in the marginal 
point s = 1 m. 

By means of the class of networks introduced in this 
work, a different type of control of the diffusion exponent 
can be realized compared to previous models. The con- 
trol parameter is not the index s as in the Levy flight, 
which is set to its marginal value s = 1 here, nor the 
prefactor /3 as in the marginal Benjamini-Berger model 
but the number fc and other possible parameters appear- 
ing in the construction procedure. The prefactor (3 in 
the distribution of lengths, as well as the shortest-path 
dimension and the random-walk dimension of these net- 
works are non-trivial functions of these parameters. 

We have also shown that by a sequence of fc = 1 ape- 
riodic networks, the normal diffusion limit (diw = 2) can 
be approached arbitrarily closely. According to numeri- 
cal results for fc = 1,2,3 uniform networks, the random- 
walk dimension decreases with the parameter fc and if 
k —^ oo, the latter quantity tends presumably to the bal- 
listic limit drw = 1- 

As opposed to the diffusion exponent that varies con- 
tinuously with the index s in case of the Levy flight, the 
set of possible values of drw which can be realized by 
these networks is discrete. Nevertheless, one can define 
uniform networks where the parameter k is not fixed but 
more than one values of fc are used in the construction 
procedure, e.g. k € {n, n -I- 1} with a fixed positive inte- 
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ger n. The random- walk dimension of such "mixed" net- 
works is expected to interpolate between that of "pure" 
ones constructed with a single k. 

In a sub-class of networks introduced in this paper, the 
dimensions dmin and have been calculated exactly; in 
some cases they have been estimated by numerical meth- 
ods. The possible exact calculation of these dimensions 
for other networks of this class or finding rigorous bounds 
on these quantities is the task of future research. 

APPENDIX A 

Let us consider a network built on an open one- 
dimcnsional lattice of size N. Connect the sites at the 
ends of the lattice (site 1 and N) with a directed link 
which the walker can traverse only in one direction, say 
from site 1 to site N with a unit rate. Let us denote the 
probability that the random walker can be found at site i 
in the steady-state by tt,;. These probabilities satisfy the 
set of linear equations 

W^J (tt, - TTj) = 0, l<i< N, (Al) 

j 

where the summation goes over the set of sites connected 
with site i and Wij denotes the jump rate from site i to 
site j; in our case wni — and Wij — 1 for other pairs 
of connected sites. Now, one can notice that Eq. 



lAll is analogous to Kirchhoff 's first law for electric cir- 
cuits. Namely, tt^ plays the role of the potential at site 
i and l/wij corresponds to the resistance of link 
Consequently, if the part of the network between site 1 
and N (except of the directed link) is replaced with a 
single (symmetric) link with a jump rate 1/f, where r is 
the effective resistance of the equivalent resistor network, 
the steady-state current J through the directed link re- 
mains unchanged. On the other hand, we may write for 
J in this simplified "circuit" , which consists of site 1 and 
site N, as well as the directed link and the effective sym- 
metric link connecting them: 

J = TTl = {tTn - TTl)f-\ (A2) 

If TV 3> 1, the effective resistance is also large and we 
have J = TTi « TTNf~^. If all the links were symmetric, 
the walker would be distributed uniformly on the network 
in the steady state. Although the link between site 1 and 
site N is directed and this leads to that the network is 
depleted on the side containing site 1, the steady-state 
probabilities far from the end 1 such as ttjv are still pro- 
portional to 1/A^ for large N. The expected value of the 
time t that the walker needs to make a complete tour on 
the network (from site N to the same site through the 
directed link) is related to the current as t — 1/J. Thus, 
it scales with the size of the network as i « ^jv^'^ ~ N^'^'' 
in the large N limit, from which we arrive at Eq. ([2|). 
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